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ABSTRACT 

We examine how the total X-ray luminosity correlates with stellar mass, stellar population, and redshift for a K-band 
limited sample of ~3500 galaxies at 0.5 < z < 2.0 from the NEWFIRM Medium Band Survey in the COSMOS 
field. The galaxy sample is divided into 32 different galaxy types, based on similarities between the spectral energy 
distributions. For each galaxy type, we further divide the sample into bins of redshift and stellar mass, and perform 
an X-ray stacking analysis using the Chandra COSMOS data. We find that full band X-ray luminosity is primarily 
increasing with stellar mass, and at similar mass and spectral type is higher at larger redshifts. When comparing at 
the same stellar mass, we find that the X-ray luminosity is slightly higher for younger galaxies (i.e., weaker 4000 A 
breaks), but the scatter in this relation is large. We compare the observed X-ray luminosities to those expected from 
low- and high-mass X-ray binaries (XRBs). For blue galaxies, XRBs can almost fully account for the observed 
emission, while for older galaxies with larger 4000 A breaks, active galactic nuclei (AGN) or hot gas dominate the 
measured X-ray flux. After correcting for XRBs, the X-ray luminosity is still slightly higher in younger galaxies, 
although this correlation is not significant. AGN appear to be a larger component of galaxy X-ray luminosity at 
earlier times, as the hardness ratio increases with redshift. Together with the slight increase in X-ray luminosity 
this may indicate more obscured AGNs or higher accretion rates at earlier times. 
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1. INTRODUCTION 

Supermassive black hole growth and galaxy growth have long 
been thought to be correlated, as supported by the Magorrian 

and Mbh — er relations (Magorrian et al. 1998; Ferrarese & 
Merritt 2000; Gebhardt et al. 2000). Both active galactic nuclei 

(AGN) activity and star formation peak at z ~ 2 (e.g., Madau 

et al. 1998; Merloni et al. 2004), and the processes appear to be 
correlated in both AGN host and normal galaxy samples (e.g., 
Canalizo & Stockton 2001; Kauffmann et al. 2003). However, 
the connection between star formation activity and black hole 
mass accretion in galaxies is not well-understood. 

Several studies suggest that AGN activity is needed to 
explain the cutoff at the bright end of the galaxy luminosity 
function and the halting of star formation in the most massive 
galaxies (Croton et al. 2006, and references therein). This theory 
is supported by the observation that AGN activity lags star 
formation by 10 7 — 10 8 yr (Davies et al. 2007; Schawinski et al. 
2009; Wild et al. 2010), and that many AGNs lie in the “green 
valley” (e.g., Schawinski et al. 2010; Xue et al. 2010) — the 

sparsely populated region between the red sequence and the 
blue cloud (Faber et al. 2007). 

However, Rosario et al. (2013b) find that out to z ~ 2 
AGNs are preferentially hosted by star-forming, rather than 
quiescent or quenching galaxies. Additionally, dust attenuation 
may redden the color of many AGN host galaxies (e.g., Aird et al. 
2012; Cardamone et al. 2010; Brammer et al. 2011), implying 
that star formation quenching by AGNs may not control the 
transition between blue and red galaxies. Instead, the connection 


between AGN and star formation activity may simply reflect 
the fact that both processes require a large gas supply, with 
possible delays between star formation and AGN activity caused 
by inefficient black hole fueling. For example, simulations by 
Hopkins (2012) and Hopkins et al. (2012) suggest that star 
formation peaks as a function of gas density, while torques 
regulating further inflow are relatively inefficient and limit black 
hole growth, leaving a reservoir of gas long after the gas supply 
for star formation is exhausted. 

While it is difficult to disentangle the contributions of star 
formation and AGN activity of high redshift galaxies from 
optical observations alone, X-ray imaging acts as a better tracer 
of AGN activity, with unobscured AGNs emitting heavily in 
both hard and soft bands, and obscured AGN dominated by 
hard X-ray emission (e.g., Churazov et al. 2002). Galactic 
X-ray emission also traces stellar activity, through low-mass 
X-ray binaries (LMXBs), whose primarily soft X-ray emission 
correlates with stellar mass, as well as high-mass X-ray binaries 
(HMXBs) and supernova, whose harder emission is reflective 
of star formation. Soft X-ray emission from hot gas is also 
apparent in the halos of massive ellipticals (e.g., Grimm et al. 
2003; Gilfanov 2004; Gilfanov et al. 2004; David et al. 2006). 

The X-ray galaxy connection has been systematically studied 
at low redshift, providing a rough census of X-ray luminosity 
with stellar mass, galaxy type, and redshift (e.g., Symeonidis 
et al. 2011; Lehmer et al. 2007, 2008). For example, Lehmer 
et al. (2007) find that the AGN fraction for massive early 
type galaxies increases with redshift, while the X-ray emission 
in less massive galaxies evolves very little out to z — 0.7. 
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Late-type galaxies evolve significantly in X-ray luminosity with 
redshift, and furthermore their X-ray emission shows a strong 
dependence on stellar mass (Lehmer et al. 2008). At higher 
redshifts. X-ray galaxy studies have primarily focused on direct 
X-ray detections (e.g., Silverman et al. 2009; Kocevski et al. 
2012; Mullaney et al. 2012; Santini et al. 2012) or on specific 
galaxy populations, such as Lyman break galaxies (LBGs), 
submillimeter bright galaxies (SMGs), and distant red galaxies 
(DRGs). For example, Laird et al. (2010) find that the X-ray 
emission is dominated by an AGN in only 15% of the SMGs 
out to z = 4 in the Chandra Deep Field North (see also 
Georgantopoulos et al. 2011), while for LBGs and DRGs the 
X-ray emission seems almost completely accounted for by 
star formation. Flowever, these different selection techniques 
target galaxies at different redshifts and stellar masses, and 
even combined are not representative of the full distant galaxy 
population. Hence, a complete census of the X-ray properties of 
distant galaxies is still missing. 

Using the high-quality near-infrared photometry from the 
NEWF1RM Medium Band Survey (NMBS; van Dokkum et al. 
2009; Whitaker et al. 2011) in combination with Chandra 
data, it is possible to systemically study the X-ray properties 
of an observed A'-band complete sample of distant galaxies. 
In this study, we utilize the ~3500 galaxy from the NMBS 
and 32 spectral energy distribution (SED) templates of Rriek 
et al. (2011), along with the overlapping Chandra COSMOS 
(C-COSMOS) field, to better characterize the star formation- 
AGN connection out to z ~ 2. With accurate Wn a and /X4000) 
measurements from Kriek et al. (201 1), indicating instantaneous 
star formation rates and stellar population ages, along with 
separate hard and soft band X-ray luminosities, we attempt to 
correlate AGN activity with star formation properties, stellar 
mass, and redshift. 

2. DATA AND METHODS 

In this work, we make use of the photometric catalogs 
from the NMBS in the COSMOS field (van Dokkum et al. 
2009; Whitaker et al. 201 1). The NMBS uses five custom NIR 
filters (three /-band and two //-band), covers 0.2 deg 2 of the 
COSMOS field (Scoville et al. 2007), and is complemented by 
FUV-MIR data sets (e.g., Sanders et al. 2007; Capak et al. 2007; 
Barmby et al. 2008; Erben et al. 2009; Hildebrandt et al. 2009). 
The photometric redshifts and stellar population properties are 
derived using EAZY (Brammer et al. 2008) and FAST (Kriek 
et al. 2009). 

We utilize the NMBS sample of ~3500 galaxies at 0.5 < 
Z < 2.0 of Kriek et al. (2011; with mass-redshift distribution 
displayed in Figure 1), and employ the spectral classification 
therein. In summary, the galaxies were divided into 32 subsam- 
ples based on similarities in their SEDs. First, the rest-frame 
photometric SEDs of all individual galaxies were compared, and 
based on the similarities in their full SED shape, it was decided 
whether or not two galaxies were analogs. Next, the galaxy with 
most analogs was identified, which together formed the first 
subsample. After the removal of all galaxies of the first subsam- 
ple, the galaxy with the next most analogs was identified. This 
procedure was repeated until nearly all galaxies were divided. 
Eventually, some galaxies were re-assigned, if they would have 
been a better match with a later iteration. Thus, our classification 
method is model-independent. 

A composite SED is constructed for each subsample by 
combining the rest-frame photometry of all galaxies in the 
bin. Each subsample contains 22-436 galaxies for which there 
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Figure 1. Mass-redshift distribution of the NMBS galaxy sample used in this 
study, with direct X-ray detections in red, and source density in gray scale. 
X-ray exposures range from 80 to 160 ks, at a sampling of ~1" pixel -1 . 

(A color version of this figure is available in the online journal.) 

is X-ray coverage as well. The final composite SED sample 
spans almost the entire galaxy population at these redshifts. 
Each composite SED shows detailed spectral features such as 
the Ha+[N n]+[S ii] and H/i+|G m| emission lines, Balmer or 
4000 A breaks, Mg ii absorption at 2800 A, the continuum break 
at 2640 A, and the dust absorption feature at 2 1 75 A (see Figure 4 
of Kriek et al. 2011 and three examples in Figure 2). As 
the sample is K-band (and not stellar mass) selected, it will 
be slightly biased toward star-forming galaxies with lower 
M/L ratios. Furthermore, as our photometric redshift code 
only uses stellar population templates, we have a bias against 
AGN-dominated SEDs. 

Although D(4000) measures stellar population age (as stellar 
opacity increases with age and strengthens absorption lines of 
ionized metals such as Ca n H and K), rather than ongoing star 
formation activity, Kriek et al. (2011) show a clear correlation 
between D(4000), Wn a , and the specific star formation rate 
(SSFR). Galaxies with shallow 4000 A breaks have high Wn a 
and are fitted by stellar population models with high SSFRs, 
while galaxies with strong breaks have low values of Wn a and 
low SSFRs. Basic properties of each SED type are shown in 
Table 1. 

In this paper, we study the X-ray properties of the different 
SED types. We make use of the C-COSMOS survey (Elvis et al. 
2009; Puccetti et al. 2009), which has an effective exposure 
of 80-160 ks, and a sampling of 1" pixel -1 . Merged images 
of 0.5-2 keV (soft), 2-7 keV (hard), and 0.5-7 keV (full) are 
publicly available, along with the point source catalog of Elvis 
et al. (2009). Each field of the C-COSMOS survey is observed 
with up to six overlapping pointings, ensuring relatively uniform 
sensitivity, despite differing point-spread functions (Puccetti 
et al. 2009). 

We obtain 60" x 60" C-COSMOS cutouts of each galaxy field 
from Kriek et al. (201 1). As a vast majority of the galaxies are 
not detected in the X-ray, we opt to stack by SED type, and 


2 


The Astrophysical Journal, 783:25 ( 8 pp), 2014 March 1 


Jones et al. 



Figure 2. Left: three examples of composite SEDs (from Figure 4 of Kriek et al. 
2011), with dashed line labeling Ha, ordered by increasing Ha strength. Right: 
corresponding 30" x 30" X-ray stacks. Green circles represent 5" apertures 
around the location of the NIR detection. External sources in this stack are not 
masked. 

(A color version of this figure is available in the online journal.) 

further subdivide by redshift and/or stellar mass. We note that 
direct detections are included in the stacks and thus the further 
analysis. We select source counts within a 5" radius of the NIR 
detection, and make background estimates by taking 30 random 
5" radius apertures in each frame, excluding 5" apertures around 
the optical source or known X-ray sources (Elvis et al. 2009). 

We then de-redshift the source according to its photometric 
redshift (or spectroscopic redshift when available), using a 
photon index of T = 1 . 1 , to determine a source luminosity, 
via the equation 

Le x -e 2 = ^xd 2 L f El -E 2 (\ + z) r ", 

where fE 1 -E 2 is the observed frame emission in the E\ — Ei 
bandpass, and cli is the luminosity distance. We explore alternate 
photon indices, as the steepness of galaxy X-ray spectra varies 
according to the strength of the various types of X-ray sources, 
and find that the observed trends in X-ray luminosity do not 
change significantly with differing T (see Section 3.2 for further 
discussion). We then determine the average luminosity of all 
sources per bin. 

To determine the error in each stack, we bootstrap resample 
the sources, subtracting a random background aperture (of the 
selected 30) from each frame. The error is the standard deviation 
of the average stack luminosities of the bootstrap simulations. 


Table 1 

Stellar Population Properties and X-Ray Luminosities 


SED Type 

A 3 

£>(4000) 

WHa log(SSFR) 
(A) (yf 1 ) 

Av log L x b log M/lo f„ M0 c 
(mag) (ergs *) (ergs *) 

i 

121 

1.95 

14.4 

-11.79 

0.4 

41 4 + 0-1 
41 -^- 0.2 

40 - 4 -°(I 2 

2 

291 

1.94 

14.4 

-11.79 

0.2 

41.6!° 0 2 3 

40.6!°j 2 3 

3 

436 

1.81 

6.2 

-24.93 

0.5 

40.5/ q 2 3 

oq 7 + 0.2 

Dy -'-0.3 

4 

211 

1.77 

2.1 

-18.35 

0.3 

41-4/oi 

40.7/° 0 2 2 

5 

95 

1.70 

24.6 

-10.59 

0.7 

41-8/.0.2 

41-0*0.2 

6 

22 

1.67 

19.1 

-11.67 

0.1 

41 7 +0 - 2 
0.3 

40 - 8 -0 2 3 

7 

72 

1.66 

20.5 

-11.29 

0.2 

41 Q+0.1 

41 O +01 
^ a - u - 0.2 

8 

64 

1.66 

38.7 

-10.49 

1.2 

41 A+0-1 
^ A -°- 0.2 

40.8/ q 2 

9 

94 

1.65 

42.7 

-10.59 

1.1 

41 3 +0,2 
0.3 

40 5 +0 - 2 
* u - d -0.3 

10 

60 

1.62 

33.1 

-10.58 

0.5 

42 2 +0,2 
^ z - z - 0.2 

41 4 + 0.2 
^ a -^- 0.2 

11 

95 

1.58 

39.1 

-10.4 

0.7 

41 Q+0-1 
4 Ly - 0.2 

41 l +01 
^ 1 • — 0.2 

12 

46 

1.58 

17.7 

-10.58 

1.4 

41 q+0-2 

41 0 +0-2 
^ 1 • — 0.3 

13 

114 

1.45 

63.3 

-10.71 

2.4 

41 8 +0 - 2 
^ 1 -°-0.3 

41 0 +0 - 2 
^ 1 • — 0.3 

14 

71 

1.45 

69.4 

-9.66 

2.6 

41 5 +01 
41.J_o2 

40 - 8 - 0 1 2 

15 

64 

1.40 

61.7 

-10.34 

2.8 

41 4+0-2 
^ a -^- 0.2 

40-7-o 2 2 

16 

155 

1.40 

67.4 

-9.92 

2.1 

41. 5^ 

40.9Z.q2 

17 

103 

1.40 

66.2 

-9.66 

1.9 

41-8-0.2 

40.2Z.Q2 

18 

95 

1.38 

86.0 

-9.34 

2.3 

41 6 +01 

"+ 1 . 0— 02 

41 0 +0-1 
^ A ' U - 0.2 

19 

33 

1.37 

86.1 

-9.66 

2.6 

41 6 +01 
^ A • — 0.2 

40 - 9 Z°o . 1 2 

20 

34 

1.37 

108.0 

-8.81 

2.6 

41 8 +03 
^ A 0.8 

41 1+0-3 
4 A - A — 0.8 

21 

139 

1.35 

79.5 

-9.15 

1.9 

41 6 +01 
0.2 

41 l +0 ‘ 1 
^ A ' A — 0.2 

22 

82 

1.34 

66.9 

-10.24 

1.7 

41 5 +0 - 2 
^ A • — 0.4 

41 l +0 - 2 
^ A ' A — 0.4 

23 

51 

1.33 

92.4 

-8.69 

2.2 

41 2 +0,3 
^ A ' — 0.9 

40 7+0 - 3 

0.9 

24 

26 

1.30 

130.7 

-9.25 

2.1 

42 1+0-2 
^•1—0.4 

41 5 +0 ' 2 
^ A - D — 0.4 

25 

159 

1.30 

94.1 

-10.71 

1.4 

41 1 +0,2 
^ 1 ’ 1 —0.4 

40.8Z°q 2 4 

26 

80 

1.30 

107.1 

-8.81 

1.9 

41 - 5/.03 

41 l +0 - 1 
^ 1 ■ A — 0.3 

27 

108 

1.27 

105.8 

-8.78 

1.6 

41-8/oi 

41 6 +0 - 2 
4 A -°- 0.2 

28 

144 

1.25 

106.2 

-9.92 

1.2 

41 8 +0 - 2 
0.3 

41 7 +0 - 2 
^ A -'— 0.3 

29 

64 

1.22 

102.8 

-8.61 

1.4 

41 7 +01 
^ 1 ’ — 0.2 

41 6 +ai 
^ a '°- 0.2 

30 

142 

1.21 

104.0 

-8.65 

1.2 

41 7+0 - 2 
^ 1 • —0.4 

41 7+0.3 
^ A ■ —0.4 

31 

101 

1.18 

109.7 

-8.65 

1.0 

41 6 +0-3 
^ 1 • — 0.7 

41 7+0.3 
^ A ■ —0.7 

32 

50 

1.13 

138.5 

-8.51 

0.9 

42 o+O - 3 
^ z - u -i.o 

42 l +0 - 3 
^ z - A -i.o 


Notes. The 4000A break (Z>(4000)) strength. Her equivalent width (WHcr), 
specific SFRs, dust attenuation (A„) are adopted from Kriek et al. The specific 
SFR (SSFR) and Ays are derived by comparing the composited SEDs to 
the Bruzual & Chariot (2003) stellar population synthesis models using a 
Calzetti et al. (2000) law, a Salpeter (1955) initial mass function, and a delayed 
exponential star formation history. The 4000 A breaks and Htr equivalent widths 
are directly measured from the composite SEDs. 
a Number of galaxies per spectral type. 
b Average full band X-ray luminosity of the stack. 
c Average full band X-ray luminosity weighted by mass of the stack. 

Stacked detections of less than 3a (as determined by the 
bootstrapped background estimates) were deemed insignificant. 
Sample stacks for three SED types of full-band X-ray data are 
shown in Figure 2. 

3. ANALYSIS 

In this section, we explore X-ray properties of the galaxy 
population at 0.5 < z < 2.0. To determine the dependence of 
X-ray luminosity on galaxy type, we seek to disentangle the 
effects of redshift and stellar mass. We attempt to remove the 
effects of mass and redshift by splitting galaxies into bins, 
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Figure 3. Average full band (0.5-7 keV) X-ray luminosity erg s vs. mass for each composite SED mass bin, color-coded by Z)(4000) and point size indicating the 
number of galaxies per bin, with la limits (left). The white circle represents galaxies not fit by any of the 32 SED templates. We additionally plot all composite SEDs 
binned by mass, rather than spectral type (right). Gray points indicate the entire sample, while colored points show the sample subdivided into three redshift bins. 
There appears to be an approximately linear correlation, with a slope of 1 (dashed line), with higher luminosities observed at higher masses and redshifts. 

(A color version of this figure is available in the online journal.) 


and eventually normalizing by stellar mass (Section 3.1). In 
Section 3.2, we look at X-ray hardness as well as expected 
X-ray binary (XRB) contribution to galaxy luminosity in an 
attempt to distinguish the sources of X-ray flux. 

3.1. X-Ray Correlations 

As each SED type probes a different mass distribution, we 
first explore the X-ray luminosity as a function of stellar mass 
by splitting each subsample into three mass bins (log M/Mq < 
10, 10 < log M/Mq <11, and log M/M 0 > 11) in the left 
panel of Figure 3. The X-ray luminosity is correlated with stellar 
mass, such that more massive galaxies emit at higher rates 
regardless of stellar age. We determine a least squares fit to 
the data via resampling, assuming a normal distribution with 
standard deviation as calculated by the bootstrap errors for each 
data point. We calculate the line of best fit for each resampling, 
weighing individual points by the number of galaxies in each 
stack, and determine the average and standard deviations of the 
fits to be logL t ~ l.l(±0.1)log M + 30(±1). 

To further examine this relation we stack all galaxies in 
10 bins of stellar mass, and show these as the large gray 
filled circles in the right panel of right panel of Figure 3, 
with errors calculated by resampling the galaxies in each bin. 
We further divide these 10 mass bins into 3 redshift bins 
(colored), and find that higher redshift galaxies have higher 
X-ray luminosities at the same mass, when compared to lower 
redshift galaxies. At z < 1 there is a clear correlation be- 
tween the X-ray luminosity and stellar mass, which is broadly 
consistent with a one-to-one relation. The X-ray luminosity 
is increasing with redshift, and the trend with mass seems 
to flatten. As we probe different mass ranges at different 
redshifts, the increase of X-ray luminosity with redshift fur- 
ther steepens the correlation between X-ray luminosity and 
stellar mass. Thus, to further examine the trend with red- 
shift and other galaxy properties, we hereafter normalize the 
X-ray luminosity by the stellar mass. 

As L x /M is roughly constant, we may estimate the typical 
accretion rate of the galaxy sample. Assuming a 1O 7 M 0 black 
hole in a galaxy of stellar mass 10 10 M o , and a bolometric 
correction of LboI ~ 15.8 L v (2-10keV) (Ho 2009), we find 
the average Eddington ratio of our sample is /-Boi/THdd ~ 10~~ 3 . 



Figure 4. Average full band (0.5-7 keV) X-ray luminosity erg s 1 vs. z for 
each composite SED bin, normalized by mass, and color-coded by D(4000) 
(see Figure 3). Points are drawn by bin size, with lcr limits. There appears to 
be at best a slight correlation between luminosity and redshift once the effects 
of mass evolution are removed. 

(A color version of this figure is available in the online journal.) 


We will argue below that accretion is the dominant source of 
emission. 

We further split the sample by redshift (z < 1.0, 1.0 < z < 
1.5, and z > 1.5; for nine total bins) and restack using the 
method described in Section 2, normalizing by stellar mass 
(such that we measure total! T/j/totaKstellar mass). Figure 4 
shows log L v (full)/M versus z for all bins combined. We find 
logL x /(M/10 10 M o ) ~ 0.8(±0.2)z + 40.0(±0.2), suggesting 
that X-ray luminosity is higher at earlier times for galaxies 
of similar mass. This trend does not change upon making a 
minimum mass cut at log M/M 0 ~ 10.6, the minimum mass 
detectable at z ~ 2. 
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Figure 5. Full band (0.5-7 keV) X-ray luminosity erg s -1 vs. Z)(4000), 
normalized by mass and corrected to the average redshift of the sample using 
the best fit L x — z relation, and color-coded by Av (points drawn by bin size, 
with la limits). X-ray luminosity per unit mass appears to be anti-correlated 
with Z)(4000), with a difference of approximately a dex between low and high 
values of Z)(4000). This may suggest higher black hole accretion rates in blue 
galaxies, or significant X-ray contribution from star formation. 

(A color version of this figure is available in the online journal.) 



Figure 6. Rest-frame U—V vs. V— J color for each SED type, color-coded by 
full band (0.5-7 keV) normalized X-ray luminosity erg s , corrected to the 
average redshift of the sample (z = 0.95) using the best fit L x — z relation, 
and scaled by bin size. The full NMBS sample with X-band S/N > 25 is 
shown in gray. Passive red galaxies (top left: bounded by the Williams et al. 
(2009) division separating quiescent and star-forming galaxies) have an average 
luminosity per stellar mass almost a dex lower than that of blue galaxies. 

(A color version of this figure is available in the online journal.) 


To remove the effects of any redshift correlation, we normal- 
ize by the above luminosity-redshift relation (as if all galaxies 
were at the average redshift, z = 0.95) after recombining all 
mass/redshift bins, and divide total stack luminosity by to- 
tal stellar mass contained in each stack, to explore trends in 
£>(4000) independent of galaxy mass and redshift. Figure 5 
shows mass and redshift-normalized X-ray luminosity as a 
function of stellar age (as measured by the strength of the 
4000 A break), with the colors indicating A v . We find the fol- 
lowing correlation between the normalized X-ray luminosity as 
a function of stellar population age: log £ X /(M/(10 10 M Q )) ~ 
— 1.8(±0.4)£)(4000) + 43.5(±0.6), suggesting that there is a 
moderate correlation between £>(4000) and galaxy X-ray lumi- 
nosity. This correlation may have caused the flattening or in- 
crease of the X-ray luminosity for the lowest mass bins for each 
redshift interval in Figure 3 (right), as these bins are dominated 
by star-forming galaxies with slightly higher X-ray luminosi- 
ties. It may also explain the flattening at the highest masses, 
as quiescent galaxies with slightly lower X-ray luminosities are 
dominating these bins. 

We also explore the X-ray properties of the 32 different 
galaxy types of Figure 5 in U—V and V—J color space in 
Figure 6, in comparison to the color of all detected NMBS 
galaxies with S/N > 25 (gray scale). This diagram is a powerful 
diagnostic in distinguishing the blue and red galaxy populations, 
and separates the red passively evolving galaxies from the dusty 
red star-formers; quiescent galaxies are primarily located in 
the upper left quadrant of this diagram (Labbe et al. 2005; 
Wuyts et al. 2007; Williams et al. 2009). The X-ray luminosity 
decreases along the star-forming sequence, and the quiescent 
galaxies have an average X-ray luminosity per unit mass almost 
a dex lower than that of the bluest galaxies in the sample. 


To first order, we observe a linear trend in galaxy X-ray 
luminosity as a function of stellar mass (Figure 3). After 
normalizing by stellar mass, we observe underlying redshift 
and stellar population age relations, finding that galaxy X-ray 
luminosity increases at higher redshifts (Figure 4), and also is 
lower for older stellar populations (Figure 5). 

3.2. Origin of the X-Ray Flux 

In this section, we attempt to distinguish between the stel- 
lar, AGN, and gaseous components of the observed X-ray 
flux. We start by comparing the measured X-ray flux with 
the expected X-ray flux from XRBs. To determine the ex- 
pected X-ray flux from XRBs, we utilize the star forma- 
tion rates and stellar masses derived in Kriek et al. (2011), 
and X-ray luminosity relations for HMXBs and LMXBs. 
Luminosity, stellar mass, and star formation rate relations 
are taken from Grimm et al. (2003) and Gilfanov (2004), 
such that £ 2 — iokeV = 6.7 x 10 39 ergs _1 SFR(M Q yr -1 ) and 
£ x (> 10 37 erg s -1 ) = 8.0 x 10 39 erg s _1 /(M st eiiar/10 11 M Q ) (see 
Gilfanov et al. 2004 and Persic & Rephaeli 2007 for further 
discussion of these relations). We plot the observed X-ray lumi- 
nosity versus expected X-ray luminosity from XRBs in the left 
panel of Figure 7. The blue galaxies are well-described by the 
X-ray luminosity of XRBs alone. However, for galaxies with 
strong 4000 A breaks, the contribution from XRBs cannot ac- 
count for the full X-ray flux. This implies that a large fraction of 
their X-ray luminosities comes from AGN and/or hot gas contri- 
butions. In the middle panel of Figure 7 we show the relative con- 
tribution of X-ray flux from non-stellar sources. This plot illus- 
trates that the non-stellar X-ray sources contribute a significant 
fraction of the X-ray flux, in particular for galaxies with larger 
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Figure 7. Left: average full band (0.5-7 keV) X-ray luminosity erg s -1 vs. expected X-ray binary (XRB) luminosity. XRB luminosity is derived from NIR star 
formation rate/mass estimates of Kriek et al. (201 1) for each composite SED bin, color-coded by £>(4000). Center: The fraction of the luminosity not accounted for 
by XRBs vs. £>(4000), color-coded by A„. Right: luminosity in excess of that predicted by XRBs, normalized to 10 10 M 0 . 1 a limits are labeled with arrows, and 
galaxies whose luminosity is over-predicted by the Ljt-XRB relations are shown at an excess luminosity fraction of 0. Luminosity in excess of that predicted by XRBs 
appears to be anti-correlated with D(4000). 

(A color version of this figure is available in the online journal.) 


4000 A breaks (luminosities of several galaxies with smaller 
4000 A breaks are upper limits). In the right panel of Figure 7 
we show the measured X-ray flux corrected for XRBs, normal- 
ized by stellar mass, as a function of 77(4000). Interestingly, we 

find that log LJM 1. 5(±0. 1)77(4000) + 43. 1(±0.2). This 

suggests that the X-ray luminosity difference between red and 
blue galaxies is almost entirely driven by non-stellar sources. 
We note that galaxies with higher values of A v may be domi- 
nated by emission from star formation (middle and right panels); 
however, these measurements are not well-constrained. 

Next we explore the hardness ratio, given by (77 — S)/(H + S ), 
where 77 is the number of hard band (2-8 keV) counts and S 
is the number of soft band (0.5-2 keV) counts (Figure 8), as a 
function of 7.7( 4000) and redshift to distinguish probable X-ray 
sources. We plot lines to determine the hardness ratio that would 
be observed in a galaxy with the listed obscuring column den- 
sities, assuming a photon index of 1.1 (solid) and 2.1 (dashed). 
We similarly plot the photon index as a function of redshift by 
combining the galaxies into 77(4000) and z bins, and measuring 
the hard and soft counts, assuming a constant hardness ratio 
(Figure 9). We find that the obscuring column density increases 
as a function of redshift, and may be higher for galaxies with 
intermediate values of 77(4000) (green). Interestingly, these are 
galaxies with high values of A v . This suggests a greater AGN 
contribution to the average galaxy X-ray luminosity at higher 
redshifts, consistent with the observed increase in total X-ray 
luminosity at the same redshifts. 


4. DISCUSSION 

After correcting for an increase in X-ray luminosity with 
both stellar mass and redshift, we find a difference of 1 dex 
between the X-ray luminosities of the highest and lowest 
SFR galaxies. To separate the effects of black hole accretion 
and stellar luminosity components, we subtract the expected 
XRB contribution and normalize by stellar mass, finding that a 
significant fraction of the X-ray luminosity is from non-stellar 
sources. This excess is highest for passive red galaxies. We note 
that both hot gas and black hole accretion may contribute to 
differences in L x . 



Figure 8. Hardness ratios vs. redshift, color-coded by £>(4000), with 
£>(4000) < 1.4 in blue, 1.7 > £>(4000) > 1.4 in green, and £>(4000) > 1.7 
in red. The sample is divided into six redshift bins, with points corresponding 
to bin size, and la upper and lower limits. Lines indicate the column density 
of the obscuring medium that would produce the given hardness ratios, where 
solid indicates T = 1.1 and dashed V = 2.1 for obscuring column densities of 
log N = 20, 21, 22, 23 (cm -2 ) (bottom-top). While there is no clear correlation 
with hardness and D(4000), the galaxies appear to be more obscured at higher 
redshift. 

(A color version of this figure is available in the online journal.) 


While hot gas likely remains an important part of the X-ray 
luminosity function out to at least z = 1.4 (Tremmel et al. 20 1 3), 
we note that its contribution to the rest- frame soft band decreases 
significantly at higher redshifts, as galactic hot gas primarily 
emits at kT < 1.0 keV. Additionally, X-ray emission from 
hot gas is unlikely directly correlated with star formation rate 
for galaxies of similar masses. The X-ray luminosities of local 
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Figure 9. Photon index as a function of redshift, as measured in 3 different 
D(4000) bins color-coded by Z)(4000), with D(4000) < 1.4 in blue, 1.7 > 
Z)(4000) > 1.4 in green, and Z)(4000) > 1.7 in red. Photon index appeal's to 
decrease as redshift increases, suggesting galaxies are more obscured at higher 
redshifts. 

(A color version of this figure is available in the online journal.) 


galaxies can vary by orders of magnitude for galaxies of similar 
stellar or black hole masses (Fabbiano 1989; Pellegrini & Ciotti 
1998; Sarazin et al. 2001; Pellegrini 2010). The cosmological 
hot gas accretion rate onto halos is predicted to plateau at z < 1 , 
however the used models are still highly uncertain (van de Voort 
et al. 2011a, 2011b). 

We thus attribute the trends in X-ray luminosity to differences 
in AGN accretion rates. Previous studies have found similar 
trends in galaxy mass and accretion rates, though we probe 
lower black hole accretion rates, averaging Froi/ThcM ~ 10 3 
(assuming the Mbh — er relation; Magorrian et al. 1998). Both 
Bongiorno et al. (2012) and Aird et al. (2012) find an increase 
in AGN luminosity with stellar mass out to z = 3 and z = 1 , 
respectively, for objects several orders of magnitude brighter 
than those in our sample. Aird et al. (2012) finds that the 
distribution of Eddington ratios is independent of stellar mass, 
suggesting that the physical mechanism responsible for fueling 
AGN is similar at all scales and times. Our results support this 
picture, with no increase in X-ray luminosity per unit stellar 
mass from 9.4 < log M/M Q < 11.6. 

Making the plausible assumption that the non-stellar X-ray 
emission is dominated by AGN, our study suggests that black 
holes have higher accretion rates in star-forming galaxies than 
quiescent galaxies, though there is significant scatter in the 
relation. This result may seem in contrast to many studies that 
found a lack of SFR-AGN connection up to z ~ 3 (Symeonidis 
et al. 2011; Cano-Diaz et al. 2012; Farrah et al. 2012; Rosario 
et al. 2012, 2013a). We note that these studies have not corrected 
for the mass evolution of their sample, and thus probe different 
galaxy populations at different redshifts. Mullaney et al. (2012) 
find that AGN host galaxies have SSFRs ~20% lower relative 
to a normal galaxy sample out to z ~ 3, but do not probe 
the lowest X-ray luminosity regimes, selecting only sources 


with L x > 10 42 . This paper explores a different regime than 
the high-SFR, high accretion rates of the bright AGNs that 
are likely associated with mergers (e.g., Santini et al. 2012). 
Low luminosity AGNs are most likely fueled by minor mergers 
and disk instabilities, and may accrete at low rates for periods 
which long outlast any intense periods of induced star formation 
(Schawinski et al. 201 1 and references therein). 

While this study probes a large range of galaxy SED types, 
masses, and redshifts, we note several limitations of this 
study. First, by stacking galaxies we are observing the average 
detection within the stack; unusually luminous galaxies could 
dominate the stacked detection. We compare to the C-COSMOS 
point source catalog (Puccetti et al. 2009), and find that in 
our stacked galaxy populations, with the exception of a few 
high-redshift bins in which one lone X-ray source is detected, 
non-detections typically provide >85% of the stack luminosity. 
This suggests that luminosities of individual bright AGNs 
do not dominate the stacks. This is perhaps unsurprising, as 
this study automatically eliminates the more unusual SEDs 
of Type I AGNs by fitting to galaxy spectra when deriving 
photometric redshifts (a vast majority of low-luminosity AGNs 
appear obscured, while those at higher luminosity are primarily 
unobscured; see, e.g., Merloni et al. 2004; Hasinger 2008; 
Treister et al. 2010). Bias from the removal of Type I AGNs 
may be most prominent for the high mass galaxy sample, as 
these galaxies are more frequently AGN hosts and have a higher 
unobscured fraction (e.g., Georgakakis et al. 2011; Fanidakis 
et al. 2012; Hopkins et al. 2012). Galaxies with SEDs that do not 
resemble any galaxy SED template (17% of the sample) average 
a full-band luminosity of 10 420 erg s -1 , significantly higher 
than most unnormalized SED types, suggesting that AGNs may 
contribute significantly to the luminosities and SED shapes of 
these galaxies. 

5. SUMMARY 

We explore the X-ray properties of ~3500 galaxies at 0.5 < 
Z < 2.0 by stacking the galaxies as a function of SED 
type (categorized by 32 SED templates in Kriek et al. 2011), 
while controlling for mass and redshift. We find that Lx 
is roughly linearly correlated with galaxy mass, and also 
increases slightly with increasing redshift. After normalizing 
by galaxy mass and correcting for the redshift correlation. 
X-ray luminosity decreases with an increase in stellar population 
ages (as measured by D(4000)). 

We utilize local XRB luminosity functions to determine the 
stellar contribution to the galaxy X-ray luminosities. We find 
that there is a significant luminosity excess for most galaxies, 
particularly those with high values of D(4000), suggesting large 
X-ray contributions from hot gas and/or black hole accretion. 

Consistent with other studies, we find that low luminosity 
AGNs show at best a mild preference for galaxies with shallower 
4000 A breaks and higher SSFRs (e.g., Xue et al. 2010; 
Georgakakis et al. 2011; Rosario et al. 2013a; Hainline et al. 
2012; Kocevski et al. 2012). First, there is a large scatter in 
X-ray luminosity from non-stellar sources for similar galaxy 
types. Second, there is no strong correlation between the X-ray 
emission from non-stellar sources and spectral type. However, 
we observe that obscuration (as determined by hardness ratio) 
appears to increase with increasing redshift as the photon index 
decreases. This is consistent with the slight increase in X-ray 
luminosity per unit stellar mass with redshift, suggesting an 
increase in AGN accretion rates in these galaxies. We note that 
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we do not probe the highest black hole accretion rates, as our 
selection methods bias against galaxies for which the continuum 
emission is not dominated by stellar light. 

While it is difficult to differentiate sources of X-ray emission 
from C-COSMOS data alone, a similar analysis of COSMOS 
radio data (such as that undertaken by Carilli et al. (2008) of 
LBGs in the field), mid-infrared data, and 6-79 keV X-ray data 
from NuSTAR could provide more direct information about the 
obscured AGN activity in these galaxies, helping to differentiate 
them from heavily star-forming galaxies that are X-ray hard. In 
addition, NIR spectra of z ~ 1.5 galaxies can aid in more 
accurate measurements of AGN demographics using emission 
line ratios [Nul/Hcr and [0 1 1 1 1 / H /! . An analogous analysis 
using the deeper Chandra Deep Field South data could also 
be performed, helping to constrain the amount of variability 
within individual galaxy SED types. 

We acknowledge the NMBS and the Chandra COSMOS 
teams for their observations and creation of the catalogs. We 
thank Ryan Hickox and Jane Rigby for useful discussions. 
T.J. acknowledges support of the National Science Foundation 
Graduate Research Fellowship. 
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